Numerical Study of Transition to Stripe Phases in High Temperature 
Superconductors Under a Strong Magnetic Field 
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The nature of spin density wave (SDW) and charge density wave (CDW) in the mixed state of 
high T c superconductors (HTS) is investigated by using the self-consistent Bogoliubov-de Gennes 
equations and an effective model Hamiltonian with competing SDW and d-wave superconductivity 
interactions. We show that there exists a critical onsite Coulomb interaction U c - For optimally 
doped sample, a two dimensional SDW and CDW modulations are induced for U <U C while SDW 
and CDW orders become antiferromagnetic (AF) and charge stripes for U > U c . These stripe orders 
are stabilized and enhanced near the vortex cores. The wavelengths of the AF stripes and charge 
stripes are found respectively to be 8a and 4o with a as the lattice constant. We show that our 
results could be applied to understand several recent experiments on HTS. 
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PACS numbers: 74.20.-z, 74.25. Jb, 74.60.Ec 



Intensive efforts have been focused on the searching 
and the understanding of spin-density wave (SDW) and 
other phases in the mixed state of high T c supercon- 
ductors (HTS) for the past several years. Experiments 
from the neutron scattering 1 , ^, H , scanning tunnel- 
ing microscope(STM) [|, £], and nuclear magnetic 
resonance(NMR) JtJ provided vital information on these 
topics. For example, according to the neutron scatter- 
ing experiment by Lake et al. a remarkable antifer- 
romagnetism or SDW appears in the optimally doped 
La 2 _ a; Sr a .Cu04 when a strong magnetic field is applied. 
Recently, Hoffman et al. || studied the LDOS in the 
mixed states of optimally doped BiaS^CaCu^Og+a using 
STM measurements, and they found that associated with 
the SDW, anisotropic charge density wave (CDW) exists 
both inside and outside the vortex cores. The coexistence 
of d-wave superconductivity (DSC), SDW and CDW or- 
ders in terms of the stripe phases was theoretically stud- 
ied in the absence of a mag netic field |, |, [l0| 



Although the competition between SDW and DSC in 
a magnetic field was previously examined jij, 13, 14|, [l5|, 
the nature of the SDW and CDW and their spatial vari- 
ations have not been addressed in such detail as to com- 
pare with the experiments. In this paper, we shall adopt 
the method described in previous papers Jll| |l7j to ex- 
amine the possible existence of SDW and accompanying 
CDW orders in the mixed state of HTS, and their nature 
in optimally doped samples. In order to simplify the 
numerical calculation, we shall assume a square vortex 
lattice for the mixed state and a strong magnetic field B 
such that A>!)>( with A as the London penetration 
depth, £ the coherence length and b the vortex lattice 
constant. Under this condition, the applied magnetic 
field B can be regarded as a constant throughout the 
sample. Our calculation is based upon a model Hamilto- 
nian with competing DSC and SDW orders and realistic 
band structure parameters. For an optimally doped sam- 
ple (x = 0.15) under a strong magnetic field B, we find 
that DSC, SDW and CDW stripe phases could be in ex- 
istence and they are pinned and enhanced by the vortex 



lattice. Numerical calculation based upon a magnetic 
unit cell of 48 x 24 lattice sites shows the wavelength of 
the SDW stripe to be 8a and the accompanying CDW 
stripe to be 4a. A phenomenological model will be used 
to explain the anisotropic SDW and CDW observed by 
experiments |2| ||. 

Let us begin with an effective mean field model in 
which interactions describing both DSC and antiferro- 
magnetic (AF) orders in a two-dimensional square lattice 
are considered. The effective one band Hamiltonian can 
be written as: 

H = ^-tijC^Cjv + ^2(Unte - (ijc^Civ 

ij.o- i.o- 

+ £(A iJ c i t T cj i + /i.c.) . (1) 



where c\ a is the electron creation operator and fi is 
the chemical potential. In the presence of magnetic 
field B, the hopping integral can be expressed as ii.j = 
t®- } exp[i-^ f*' A(r) • dr] where t?j = t for the nearest 

neighboring sites while the next-nearest neighbor 

hopping t?j = t' . The superconducting flux quanta de- 
notes as $o = h/2e. Here we choose Landau gauge 
A = (— By, 0,0) with y as the y-component of the po- 
sition vector r. The two possible orders in cuprates 
are SDW and DSC which have the following definitions 
respectively: Af DW — U{chc^ — c ij. c u) an d Aij = 
Vdsc (ciT c jl — c ii c jt)/2- bi the above expressions, U and 
Vdsc are respectively the interaction strengths for SDW 
and DSC orders. Vdsc, which gives rise to the d-wave su- 
perconductivity, may come from all possibilities including 
AF fluctuations and electron phonon interactions. The 
mean- field Hamiltonian (1) can be diagonalized by solv- 
ing the resulting Bogoliubov-de Gennes (BdG) equations 
self-consistently 
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where the single particle Hamiltonian Hf. 
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(Uriig- — /i)<5ij, and and n,-f 
E>Hl 2 (l - /(£»)). A y - 

u il M if ) tan ^ ( 2fc2*T ) ' w ^ n /(-^) as tne Fermi distribu- 
tion function and the electron density ni — + rijj. 



The DSC order parameter is defined at site i as AP = 
(AP +ex , ; + Af_ ex , - AP i+ey - AP ; _ ey )/4 where Aft = 

"(ri+rj)/2 



Aijexp{i-^ J^ r ' +r J^ z A(r) • dr] and e x y denotes the unit 
vector along (x, y) direction. The main procedure of self- 
consistent calculation is given below: For a given initial 
set of parameters and A; j, the Hamiltonian is nu- 
merically diagonalized and the electron wave functions 
obtained are used to calculate the new parameters for 
the next iteration step. The calculation is repeated un- 
til the relative difference of order parameter between two 
consecutive iteration step is less than 10~ 4 . By varying 
the chemical potential, one obtains solutions correspond- 
ing to various doping concentrations. 

In the following calculation, the length and energy are 
measured in units of the lattice constant a and the hop- 
ping integral t respectively. Here the next-nearest neigh- 
boring hopping integral is chosen to be t' ~ —0.2 to fit the 
band structure of HTS. It needs to be pointed out that 
the induction of internal magnetic field by the supercur- 
rent around the vortex core is very small as compared 
with the external magnetic field, so that the uniform 
magnetic field distribution is a valid approximation. We 
follow the standard procedures [|l6| |l^] to introduce mag- 
netic unit cells, where each unit cell accommodates two 
superconducting flux quanta. Periodic boundary condi- 
tion is imposed in calculation. The related parameters 
are chosen as the following: For optimal doping x = 0.15 
(or the electron doping rif = 0.85), the DSC coupling- 
strength is Vdsc = 1-0, the linear dimension of the unit 
cell of the vortex lattice is chosen as N x x N y — 48 x 24 
sites. 

First let us choose U — 2.2 such that AF order is com- 
pletely suppressed at zero field. Our calculation is per- 
formed at low temperature. The spatial variation of DSC 
order parameter Af is plotted on a 24 x 24 lattice in Fig. 
1(a) with the vortex core situated at the center where the 
DSC order parameter vanishes. By comparing it with the 
vortex structure of a pure DSC, the size of the vortex core 
here is noticeably to be enlarged. Fig. 1(b) displays the 
spatial variation of the induced staggered magnetization 
of SDW order as defined by M? = (-1)* Af DW /U. There 
the SDW order xists both inside and outside the vortex 
cores , and exhibits isotropic two-dimensional behavior 
with the period 8a along both x and y directions. Its 
magnitude reaches the maximum value at the vortex core 
center. The DSC and SDW orders coexist throuhout the 
whole sample. The appearance of the SDW order around 
the vortex cores strongly affects the spatial profile of the 
local electron density distribution, which can be repre- 
sented by a weak CDW as shown in Fig. 1(c). The re- 
markable enhancement of electron density (or depletion 



of the hole density) is presented at the vortex core center. 
The variation of the electron density outside the vortex 
core shows weak oscillation. Different two-dimensional 
SDW and CDW structures have also been obtained for a 



different set of band parameters |18 



Next we perform the calculation at very low tempera- 
ture for U — 2.4, the obtained results are fundamentally 
different from those for U = 2.2, and they are presented 
in Fig. 2. In Fig. 2(a) we plot the spatial variation of 
DSC order parameter. It is clear that y-axis oriented 
stripe like structures appear in Af 3 with a weakly modu- 
lation period of 4a. The size of vortex core is further en- 
larged and elongated along the y axis than in Fig. 1(a). 
The SDW order and CDW order are displayed in Fig. 
2(b) and Fig. 2(c), respectively. The SDW order be- 
haves like almost uniform AF stripes oscillating with a 
wavelength of 8a. The vortex core is always pinned at 
one of the ridges of AF stripes where the AF order is 
stronger than those at other sites. The spatial modu- 
lation of CDW order also exhibits quasi-one-dimensional 
charge stripes behavior with a wavelength 4a, exactly half 
of that of the SDW along the x direction. The above nu- 
merical results are checked by three different set of initial 
parameters n\ a and Ai j and the iteration processes have 
been carried out for more than 500 steps to achieve the 
required accuracy. The above results for finite B indi- 
cate that there exists a critical point U c ~ 2.25 between 
U = 2.2 and U = 2.4, such that isotropic two-dimensional 
SDW and CDW may be induced when U < U c and they 
become stripe like structures when U > U c . As we shall 
show below that the AF stripes and charge stripes ob- 
tained here could be very relevant to experiments per- 
formed on the optimally doped BSCCO. It is likely that 
some of the optimally doped and underdoped HTS could 
be close to this critical region. At B = we found no 
two-dimensional SDW and CDW regardless of the value 
of U. 

Before comparing with experiments, we would like to 
point out that the stripe phases oriented along x- and y- 
directions are degenerate in energy. In order to compare 
with the experiments, we shall assume that in certain 
domains of the sample the stripes are further pinned by 
some defects which makes their orientation along the y- 
direction more energetically favorable than those along 
the x-direction. Setting the energy difference between 
the x- oriented and y-oriented stripe phases to be SE > 
and defining rj = exp(—SE/T) with T as the tempera- 
ture, then the statistical probabilities for y-oriented and 
x-oriented stripes to appear at temperature T are re- 
spectively 1/(1 + 77) an d v/0- + 7 l)- The measured order 
parameter should be O — [0(y) + nO(x)]/(l + rj) with 
0(x) [0(y)] representing the order parameter for one of 
the x- (y-) oriented DSC, AF and charge stripes. The 
combined results are shown in Fig. 3 for 77 = 0.5. The 
spatial variation of the combined superconductivity order 
parameter is presented in Fig. 3(a). Here the SDW (see 
Fig. 3(b)) and CDW (see Fig. 3(c)) have anisotropic two 
dimensional structures and respectively with periodici- 
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ties 8a and 4a, in good agreements with the observations 
of Lake et al. H and Hoffman et al. ||. Our results 
also predict that as T approaches zero or rj = 0, AF and 
charge stripes should show up, and at higher T or 77 close 
to f , the observed SDW and CDW should become more 
isotropic. It is useful to point out that the way to explain 
these experiments may not be unique. As suggested by 
Howald et al. p9[ , if x-oriented and j/-oriented stripes 
are pinned in two nearest neighboring domains, the prox- 
imity effect may cause these two different oriented stripes 
permeating each other and make their spatial distribu- 
tion to look more two- dimensional like. So far we are 
not able to numerically simulate this situation. 

The features exhibited in Fig. 2 and Fig. 3 are very 
robust when U > U c . We believe that the optimally 
doped BSCCO sample should have a U > U c . While the 
STM experiments || were performed at B ~ 7 Tesla, our 
numerical results based upon the 48 x 24-sites calculation 
corresponding to a magnetic field B ~ 27 Tesla. For a 
realistic comparison with experiments, one needs to do 
a 92 x 46 -sites calculation, that is beyond our current 
computing capability. In a weaker field B ~ 7 Tesla, we 
expect that the SDW and CDW structures with periods 
8a and 4a should still exist in the neighborhood of a 
vortex core. Away form the vortex core, the situation 
should be corresponding to B = case. 

We have also done a numerical study for B = with 
U = 2.4 and x = 0.15, and our calculation indicates that 
both the stripe phase and the uniform DSC phase (no 
SDW and CDW) could show up depending on the initial 
input parameters. Since U — 2.4 is quite close to U c , 
the free energy difference between the stripe phase and 
the uniform DSC phase is estimated to be very small, 
which suggests that the experimental observed phase at 
finite temperature should come from a superposition of 
these two configurations. This consideration would dra- 
matically reduce the amplitudes of the SDW and CDW 
stripes measured by experiments. But when a magnetic 
field B is applied, the stripe phase is the only solution 
regardless of the initial parameters. This result implies 
that the SDW and CDW orders are stabilized and en- 
hanced near the vortex core within a distance of several 
coherence lengths. Away from the vortex core they are 
somewhat suppressed. Of course, the stripe phase order 
at B = could be stabilized by the presence of defects 
and it can also be strengthened by a larger U. The en- 
hancement of AF order near the vortex core in a strong 
magnetic field is consistent with neutron scattering ex- 
periments Jl|, H . Here The manifestation of stripe phases 
in optimally doped sample seems to be against common 
consensus. But our stripes are small one-dimensional 
SDW and CDW modulations in a d-wave superconduct- 
ing background. This is different from what was pro- 
posed originally by Emery and Kivelson for under- 
doped sample where the AF phases are insulators (no 
hole regions) and charge stripes are conductors (rich hole 
regions). There is no physical reason to forbid our stripe 
like modulations appearing in optimum doped HTS. The 



observation of charge stripes in very recent STM experi- 
ments Jl!| at B = seems to suggest that stripe phases 
may indeed be present in optimally doped BSCCO sam- 
ples. 

For the purpose to have a better understanding of the 
doping effect, the underdoped case (x = 0.10) is exam- 
ined. We found that the spatial distribution of SDW and 
CDW still exhibits the stripe-like behavior as shown in 
Fig. 2, and the periods of SDW and CDW now change 
respectively to 12a and 6a. We expect that the periods 
could even become 16a and 8a, when the doping level is 
further reduced. These results are in qualitative agree- 
ment with the observations of Lake et al. || where a mag- 
netic field induced AF order was observed in underdoped 
sample. A detailed study of the doping dependence will 
be reported elsewhere. 

Finally, we would like to point out that Zhu et al. 
I^rjdid almost the same calculation using slightly larger 
U = 2.5 and x — 0.16, their results yield isotropic two- 
dimensional SDW and CDW, which are very different 
from the AF and charge stripes obtained by us. We 
have carefully checked their calculations and found that 
if they increased the number of iteration steps sufficiently 
to achieve higher accuracy, their two-dimensional SDW 
and CDW checkerboard patterns would have evolved into 
our stripe like structures as shown in Fig. 2. 

In summary, the stabilization and the enhancement 
of AF and charge stripes in a d-wave superconductor 
near a vortex core are numerically studied by a mean 
field Hamiltonian. We found the wavelength of the AF 
stripes to be 8a and that of charge stripes to be 4a. 
Assuming that the degeneracy of the x-oriented and y- 
oriented stripe phases is broken by some defects, we show 
that the observed spatial variations of SDW and CDW 
are anisotropic and two-dimensional, with wavelengths in 
good agreement with experiments. We also would like to 
emphasize that our self-consistent mean field BdG equa- 
tion calculation tends to overestimate the stability of the 
static AF order. In order to partially overcome this defi- 
ciency, we choose the onsite Coulomb interaction U = 2.4 
to be somewhat smaller than that for the standard Hub- 
bard model. The effect due to the dynamic SDW has 
not been included in this study. Any attempt to do the 
present type of calculation by including this effect is quite 
difficult and has to be confined to much smaller mag- 
netic unit cell. This would make the comparison with 
experiments difficult. As to whether a static AF order 
could exist in optimally doped sample is still a subject 
for debate. The coexistence of charge stripes with DSC 
at B = observed by very recent STM experiments 
at optimal doping may indicate that static AF stripes 
are also in presence. For a pure DSC, the local density 
of states (LDOS) at the vortex center is well known to 
have a broad peak around E = Q. However, the van- 
ishing LDOS at E = near the vortex core observed by 
STM experiments [g |] for YBCO and BSCCO has been 
understood in terms of the presence of SDW in a d-wave 
superconductor pM , an indirect evidence of the existence 
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of the static AF order in optimally doped HTS samples. 
In view of all these and the favorable comparison with ex- 
periments, The qualitative feature of our results should 
still remain even in a more refined theory. 
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FIG. 1: Spatial variations of the DSC order parameter Af 
(a), staggered magnetization M? (b), and electron density rii 
(c) in a 24 x 24 lattice. The size of a magnetic unit cell is 
48 x 24, corresponding to a magnetic field H — $o/(24 x 24). 
The strength of the on-site repulsion U = 2.2 and the averaged 
electron density n — 0.85. 



FIG. 2: Spatial variations of the order parameters Af 3 (a), 
M? (b), and rii (c). The strength of the on-site repulsion 
U — 2.4. The other parameter values are the same as Fig. 1. 



FIG. 3: Spatial profiles of the combined x- and y- oriented 
stripes with the order parameters AP (a), Mi (b), and n\ (c). 
Only the 20 x 20 lattice with the vortex core at the center is 
plotted in (c). The mixing factor is chosen as r\ — 0.5. The 
other parameter values are the same as Fig. 2. 



Fig. 1(a) Chen et al. 
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Fig. 2(a) Chen et al. 




Fig. 3(a) Chen et al. 




Fig. 3(c) Chen et al. 




